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Abstract 

The oxygen 02pa and copper Cu4s and Cu3dj.2_y2 orbitals are involved in 
a simple LCAO model for determination of the conduction band and the 
oxygen-oxygen hopping is considered as a small parameter with respect to 
the transition amplitude between nearest neighbours. The traditional Cooper 
pairing is obtained by taking into account the double-electron exchange be- 
tween the nearest neighbours within the two-dimensional Cu02 plane. The 
equation for the superconducting gap is obtained as a result of the standard 
BCS treatment. It is shown that the order parameter could have either s-type 
or d-type symmetry depending on the ratio between the transition amplitudes. 
This model allows understanding the experiments reporting a vr-shift of the 
Josephson phase indicative for a d-type gap symmetry as well as the observed 
s-type in the case of strongly irradiated samples. 

PACS numbers: 74.20.Fg - BCS theory and its development; 74.72.-h - High- 
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I. INTRODUCTION 



The discovery of high-Tc superconductivity has brought significant interest into this 
field and triggered many intense investigations during the last decade. In the centre of 
them is the question about the determination of the basic parameter characterizing this 
phenomena - the superconducting order parameter. Recent experiments on angle-resolved 
photoemission spectroscopy (ARPES) gave fast increase of the quantitative results for the 
Fermi surface and for the type of the angular dependence of the order parameter as well. 
As a result of these studies, for YBa2Cu307 and Bi2Ca2SrCu208 symmetry of the type 
cospx — COS Py is often assumed. This assumption has been independently confirmed by the 
experiments on Josephson junctions for YBa2Cu307 0. On the other hand, deviation from 
the simple ci-case was observed in the experiments with strongly irradiated samples As 
suggested by Abrikosov and Pokrovsky and Pokrovsky that could be realized by the 
reduction of the (i-channel and domination of the s-part of the electron-electron interaction. 
The Nd2-a;Ce2,Cu04_5, for example, has similar symmetry of the superconducting gap. 

The purpose of this paper is to derive analytical expression for the interaction involved 
in the standard equation for the superconducting gap by successive BCS treatment of the 
double-electron exchange between the nearest neighbours (NN) and the next nearest neigh- 
bours (NNN) in the Cu02 plane. The matrix elements of the interaction, being a sum of s- 
and (i-symmetry terms, and the hmit cases leading to simple s- or d-tjpe gap are discussed 
as well. 

The eigenfunctions and eigenvalues for the conduction band found by solving of 
Schrodinger equation are used to obtain the momentum representation of the superexchange 
interaction. The successive BCS scheme applied to the latter leads to equation for the BCS 
gap. The interpretation of the exact result in the limit cases of strong hole and electron dop- 
ing is discussed in Sec. |V|. In conclusion the fitting of our results to the recent experimental 
data is considered. 

II. MODEL 

Following the ideas of quantum chemistry we shall use a tight-binding (TB) method 
to obtain the electronic band structure of layered cuprates. To this purpose we consider the 
atomic orbitals related to the Cu4s, Cu3dx2^y2, 02po- states. Denoting with R^,., Hy, Rcu the 
positions of O^, Oj^, and Cu atoms in the Cu02 plane, with oq the in-plane lattice constant 
and n - the unit cell index, the wave function within the adopted here liner combination of 
atomic orbitals (LCAO) approximation reads as 

V'lcao (i") = ^n,ai^02p^ (r - - R^) + >"n,aV'o2p, " ""O " ^y) (2.1) 

+Sn,a^cu4s (i" - nao - Rcu) + Dn,a^cu3d , , (r - nttQ - Rcu), 

where the coefficients X^^a, ^n,^, 'S'n,^, and -Dn,a are the amplitudes for the n-th unit cell. The 
building of TB Hamiltonian in the terms of second quantization is reduced to replacing these 
amplitudes by creation and annihilation operators satisfying the anticommutation relations 
of the type |Xn,Xj„i = 5n,m, |^n, ^j^^ = 0. Further introduce the notations Is for the 
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amplitude of the transition between the Cu4s and 02pa and Id - between CuSd and 02p^ 
orbitals. The Hamiltonian giving the band structure of CUO2 plane, which incorporates the 
oxygen-oxygen hopping amplitude t, has the form 



—Id {—Dn + Dx+l,y) + £2pa:Xri\ 

^ (-^n — Xx-l,y — Xx,y+l + Xx-l,y-{-l) 
—Id {Dn — -Dx'+l,j/) + £^2pj;^n] 



(2.2) 



where the single-site energies of 02^^, 02py, Cu4s and Cu3dx2_y2 orbitals are denoted by 
^2p^, £2py, ^is and £3^ respectively; the energy is measured from the oxygen 2p level, i.e. it is 
assumed £2p^ = £2py = £2p = 0. Hence, the energies of the copper orbitals are = £4^ — e2p 
and ed = e^d — £2p- Now using the Bloch waves 



Y 

Y 



_J_ y- ip-n (_ipipx/2\ V v] 

_J_ „ip-n ( _\f^ipy/2\ y yf 
^Z^K y ic ) ^p,a, -'n,o 



/N 



p,ai 



Si 



p 

^2^e ^le ; V,a' 



(2.3) 



nt = ^ p-ip' n nt 



where p is dimensionless momentum {px,Py) € (0, 27r) and taking into account the relation 
I]n6'^P~P = ^p',p the Hamiltonian, Eq. ( p.2| ), is reduced to the form 



where 



'4^p,a — 



D 



p,a 
p,a 



p,a 



p,a 



IsS^ 
Is'^Y 

-hs^ 





ed) 



(2.4) 



The notations used here are after Andersen et al. 



(4 + 4)^' ^ 



2sin(px/2), = 2sin(pj^/2), 



cospx)/2 and y 



cos Py)/ 2. 



To find the energy spectrum of the TB Hamiltonian one can employ the method described 
in Ref . ^ . Its essence comprises in extracting an effective oxygen part of the TB Hamiltonian 
by eliminating the metallic amplitudes (a procedure also known as Loewdin perturbation 
technique). In our case these are Sp and Dp which read as 
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Dr. 



{s^Xp SyYp). 



(2.5) 



Hence we obtain 2x2 matrix problem. The effective oxygen 2x2 Hamiltonian takes the 
form 



H, 



(O-O) 
eff 



Ho + Vt, 



Hn 



--Dcff I 



Vt = -t 



eff 







(2.6) 



where B^s = 47 (^s - e) + Ij/ (e^ - e) and Us = t - 2Jj/ (e^ - e). 

To solve the eigenvalue problem for H^^~'^^ we assume that tes B^s and will use 
perturbation theory with respect to the small parameter r = teffZ-Bes- In zeroth order 
approximation we have 



rl°) = 0, 



c(o)) = 1 



(2.7) 



For our purposes we will consider the first order correction with respect to |c) and Ec They 
are given by (see for example Ref. 0) 



1^ / ^(0) ^(0) 



(2.8) 



.(1) 



(p) = (c(°)|\4|c(°)) = 2t 



eff 



X Y 

+ 

X ' y 



One can readily obtain the required matrix element by using Eqs. ( p.6|) and ( p.7|) 



(6(o)|i/|c(o)) = -t 



Therefore, according Eq. ( p.8|) , the first correction to the |c) vector takes the form 



r 



Now substituting Eq. 

get 



]5[), for the conduction band in (r <^ l)-approximation we finally 



c) ~ - 

s 



4) \ 



2/sT / 2 _ „2 ^ 

es-e'^x^YK^x ^Y^ 



\ el-e'^^x^Y 



1 -(4-4) ' 



(2.9) 
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£e(p) = 4teffX— X- ^2.10) 

2x ~^ 2y 

The last two expressions are used to derive in the next section the four-fermion term which 
describes the interaction between electrons leading to attraction. 



III. THE HEITLER-LONDON INTERACTION 



In order to describe the effective interaction between electrons we shall start here from 
two-electron exchange Hamiltonian. The underlying idea of a double electron exchange has 
been considered, for example, in Refs. [p!0|-[T3[| and the original Heitler-London's considera- 
tions in the theory of H2 molecule consist in involving a double electron exchange amplitude 
that takes into account the correlated hopping between neighbouring atoms [Q. 

In the case of CUO2 plane the transitions between Cu4s, Cu3dx2_y2 and 02p„ must be 
taken into account. The 2po- 4s transition amplitude is denoted by Jsp in the following, 
and Jdp stands for the 2p^ 3d hopping, respectively. In order to complete the investigation, 
started in Refs. []T3|,|IB[, here we will not take into account the 0-0 hopping amplitude. 
Consequently, the four fermion interaction reads as 



n,a,/3 



'^'^n,l3'^x+l,y,a-^n,aDx+l,y,l3 + Yl,i3'Dx,y+l,ayn,aDx,y+l,f^ | 



t ct 



+ Jc 



dp 



(3.1) 



Each term could be compared to the corresponding one for H2 molecule in Ref. [jT6[, ifnL 
The direct substitution of the transformations below in the interaction Hamiltonian 



t 

n,/3 



i^i:e'P-(-ieW2)(fx)c^,,, Yl, 



/27V ^ e 



^Ee-'P'- (ie-'^i 
p 



1 2J,r„-iq'.n^ X^y(4-4) J 



/2N 'H e 

q 



(3.2) 



iqn 



^2N ^ ed-e 



2s^ 



-iq n 



'2N ^ ed-e 



2s ^s. 



^ X Y ' 



2s2 
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leads to the momentum representation of the interaction. Here we shall suppose e <^ e^, 
and therefore q = — e ~ e^. During the calculations we have used the equality 



E -ip' n -iq' n ip n iq n _ X , , 
e e e e — Op'+q',p+q, 

-' ^ n 

and thus we have a sum over four momenta which satisfies the quasimomentum conservation 
law. In the case of space homogeneous order parameter and currentless equilibrium state 
we must take into account only the terms with zero momentum p + q = p' + q' =0, or 
q = — p, q' = — p', and this leads to simplification of the result for ifint 

1 



int 



2N 



(3.3) 



p,p' ,a,P 



where 
^(P',P) = \ Jsp 



rV(p')a(p) 



dp 



ra (p, 



(3.4) 



X 



siip')si{p') 4(p)4(p) 



s(pO 



sip) 



^ cot{p'j2) cot(p J2) + ^ cot(p;/2) cot(p,/2) ) . 



and 



X Y 



We consider the case where the influence of the odd in p^, and Py terms is negligible. This, for 
example, holds for the conventional superconductors, described by the BCS theory and 



could take place in the layered cuprates as shown by the experiments on Bi2SrCa2Cu208 JTS 
In the next section we consider the possibilities provided by the reduced kernel 



^hl(p',p)= J. 



sp 



+ J. 



dp 



rV(p')a(p) 



X 



4(p')4(p') 4(p)<(p 



ra (p, 



(3.5) 



IV. THE BCS SCHEME 



Consider now Vhl(p') p) involved in the self-consistent BSC calculation |[T^ of the super- 
conducting gap. Following the method described in this fundamental work and notations 



from Ref. |19] we obtain the following expression for the order parameter 



A(P') 



1 

2iV 



Enp',p) 

pe£ 



E(p) = VA2(p)+r/2(p), 



(4.1) 
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where 77(p) = p^/2m — /i, with /i being the chemical potential of the electrons. The renor- 



malization procedure [19| tells us that the summation in Eq. (4J.) is only over a narrow 
energy interval along the Fermi surface contour (FS £), i.e. Ep — Tiuj^ < eip) < E-p + hu^, 
where the cut-off parameter of the sum hu^ is found to be hu^^ ~ E-p/ 2. For more details 



about the calculations see for example Ref. 
order parameter symmetry is studied as well 
for the constant energy contours (CEC) is given by e. 



20[| , where the influence of the impurities on the 
In the case of layered cuprates the equation 
p) = Ep, where £c(p) is given by 



Eq. ( p.lOp . The result for Vhl can be further simplified if we introduce the dimensionless 
Fermi energy measured in units of the conduction band width w 



^clPj 



W 



1 44 



2x^ 2y ^-^x^ 



const, 



where w = Ates is the bandwidth. Thus we get 



(4.2) 



^hl(p',p) = IQel \ Jsp 



+ Jc 



dp 



slip') -slip')) (4(p)-4(p; 



s[p' 



2-Ta\p: 



(4.3) 



V. DISCUSSION 

To gain further knowledge on the gap symmetry it is straightforward to examine Eq. ( |4.1|) 
for particular choices of the interaction parameters entering Eq. ( ^4.3| ) for which certain 
plausible limit cases occur. Thus, for instance, if the 3d amplitudes dominate the transitions 
between the Cu and O orbitals as considered in Ref. []2T|, one would have Jgp <^ Jdp and 
therefore 



where 




Thus we have a separable Hamiltonian of the interaction leading to an implicit analytical 
solution for the gap 




After separating the angular dependence in Xsip) and introducing the so called order pa- 
rameter SgiT) at finite temperature T, we have 



7 



A(p)=xs(p) S,(T), 
where, in compliance with Eq. ([4 .11) 



tanhfHE^i) (5.2) 



E(p') 



This expression has the standard BCS form for a scalar type gap [jl^. In this case, Jsp <^ Jdp, 
we have the angular dependence 

A(p) oc Xs{p) > 0, 

and therefore, as we use r <^ 1, it exhibits s-type symmetry, A ^ const. Such a possibil- 
ity exists in strongly irradiated samples (see also Refs. and references therein) or 
Nd2-xCe2;Cu04_5. 

Consider now the opposite case Jdp ^ Jsp- In the separable kernel obtained 



V^L ~ le^y^p [ f ) r\d{p')Xd{p), (5.3) 



the Xd{p) function has the form 

Xd{p) = slip) - slip) = -2 (cosp^ - cospy) , 

which yields the so called d-tjpe gap anisotropy. Now the gap equations Eqs. ( |4.1|) , ( ^.1|) 
and ( |5.2| ) take the form 

^(p) = Xdip) 2d(T) oc cos(p^) - cos{py), 

^ 16/ ( Is\^ 2( /\ t^nh{E{p')/2T) _ . 

2^iOJ^p£„ l^^j r 2^Xrf(^pj ^(p,) - i, (54) 



E{p) = Jxl{pmT) + v'{p). 



Here the angular dependence is carried by the fragment s —s, the nodes of the gap are now 
situated at the points = ivr ± py, p^ = ^Py, i-e. along the diagonals of the rounded FS 
square = ec{p)/w = const, in case of < 0.38 as it is for the hole doped YBa2Cu307 and 
Bi2Sr2CaCu208. Such a location of the nodes of A(p) is in accordance with that reported 
in Ref. |jl8| where at last the methodologically important pair of parameters {u, v) of the 



theory [|T9| have been measured by ARPES. 

To bridge the current discussion and the experiment we employ Eq. (|5.4|) to fit the recent 



experimental data by Ding et al. |TB| and the result is shown in Fig. |I]. Since A(p) lives 
only on the FS contour a simultaneous fit to both the gap and FS can be achieved by simply 
projecting the gap curve onto the {px,Py) plane. As clearly seen in Fig. |l|, the adopted here 
simple TB self-consistent model remarkably reproduces the experimentally observed A(p) 
anisotropy in Bi2Sr2CaCu2 08. 
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VI. CONCLUSIONS 



In the preceding sections we have given an account of a model in which different cases 
of superconducting gap symmetry occur upon 'passing' various sets of parameters as an 
'input'. 

Having obtained the underlying microscopic mechanism of high- Tc superconductivity 
we hope that subsequent investigations on different materials will finally determine the 
parameters entering the interaction Hamiltonian so that its structure be enough to explain 
the experimentally observed different types of superconductivity. Let us note that now in a 
decade of investigations it is not yet firmly recognized whether it is a standard BCS scheme 
or some kind of exotic interaction that gives rise to high-Tc superconductivity. 

In conclusion we stress that within the framework of the suggested model not only the 
results that give s-type symmetry is easily interpreted, but also the recent experiments 
on the Josephson vr-shift in YBa2Cu307 and the ARPES study of Bi2Sr2CaCu208 and 
YBa2Cu307 This model makes use of ideas having their origin in the quantum 

chemistry, quantum field theory and gives, by itself, a successive microscopic derivation of 
the interaction Hamiltonian ifjnt of the BCS theory. Moreover the interpolation formulae 
used to fit the experimental data for the Fermi contour and the angular dependence of the 
order parameter, for instance, are obtained as a simple result within the framework of the 
traditional band picture and the BCS scheme. We consider that it is most unlikely the 
same analytic interpolation formulae to be successively derived by an alternative theoretical 
model, i.e. model using only college trigonometry and exhibiting textbook-like behaviour. 
In this sense the theory of superconductivity repeats the development of quantum electrody- 
namics from half century ago and we could see the victory of traditionalism in the decadent 
theoretical physics at the end of the 20-th century. 
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FIG. 1. The analytical 'cos(pxO — cos{py)' fit (thick solid line), according to Eq. (5.4), to the 
experimental data by Dingh et al. [18] for the gap anisotropy in Bi2Sr2CaCu208 (the dots). The 
fitted Fermi surface line C (horizontal rounded square) is obtained by projecting the gap curve 
onto the {px,Py) plane. The notations for the high-symmetry points of the Brillouin zone are those 
standard for simple square lattice. 



11 



